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Abstract 

Given a parametric polynomial ideal /, the algorithm DISPGB, introduced by the 
author in 2002, builds up a binary tree describing a dichotomic discussion of the different 
reduced Grobner bases depending on the values of the parameters, whose set of terminal 
vertices form a Comprehensive Grobner System (CGS). It is relevant to obtain CGS's 
having further properties in order to make them more useful for the applications. In 
this paper the interest is focused on obtaining a canonical CGS. We define the objective, 
show the difficulties and formulate a natural conjecture. If the conjecture is true then 
such a canonical CGS will exist and can be computed. We also give an algorithm to 
transform our original CGS in this direction and show its utility in applications. 

Keywords: canonical discussion, comprehensive Grobner system, parametric polynomial 
system. 

MSC: 68W30, 13P10, 13F10. 

1 Introduction 

There are many authors [BeM IBeWe931 IDeM IDoSeSt061 IDu95l IFoGiTrOll lGi87l IGom02l 
IGoTrZaOOl lGoTrZa051 IBeMcKa97| lKa97l |Kap95j IMaMoOB} lMo02l IMor97} IFeM ISaSu03l 



SuSa06[ Sc9l| ISi92j IWe92[ IWe03t IWi06] who have studied the problem of specializing para- 
metric ideals into a field and determining the specialized Grobner bases. Many other au- 
thors [CoM IEm991 IGoB,e931 IGuOr041 lMo95l IMoM [R^OO] have applied some of these 
methods to solve concrete problems. In the previous paper [Mo02| we give more details of 
their contributions to the field. In the following we only refer to the papers directly related 
to the present work. 

Let I C K [a] [x] be a parametric ideal in the variables x = x%, . . . , x n and the parameters 
a = at, . . . ,a m , and y% and >~a monomial orders in variables and parameters respectively. 



*Work partially supported by the Ministerio de Ciencia y Tecnologfa under project MTM 2006-01267, 
and by the Generalitat de Catalunya under project 2005 SGR 00692 
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Denote A = K\a]. Weispfenning [We92] proved the existence of a Comprehensive Grobner 
Basis (CGB) of / and gave an algorithm for computing it. It exists a modern implementation 
in REDUCE of CGB algorithm due to T. Sturm et al. [DoScSt06]. Let K be a computable 
field (for example Q) and K' an algebraically closed extension (for example C). A CGB of 
I C A[x] wrt (with respect to) the termorder ^ x is a basis of / that specializes to a Grobner 
basis of cr-a (I) for any specialization o~a ■ K\a]\x] — > K'\x], that substitutes the parameters 
by values oo G K' m . 

In most applications of parametric ideals the related object called Comprehensive Grobner 
System (CGS) is more suitable. A CGS of the ideal / C A[x] wrt >~% is a set 

CGS(J, Hr) = {(Si, Bi) : 1 < i < s, Si C K' m , Bt C A[x], \Ji Si = K' m , 
Vao £ Si, o~a {Bi) is a Grobner basis of 0"a o (-O wrt >~x}- 

The sets Si are often called segments and it is always assumed that they are constructible 
sets. A CGB is a special CGS with a unique segment K' m . In a CGB the polynomials in 
the basis are faithful, i.e. they belong to /. Further properties are required to obtain more 
powerful CGS. 

Definition 1 (Disjoint CGS). A CGS is said to be disjoint if the sets Si form a partition 
of K' m . 

Definition 2 (Reduced basis). A subset B C A\x] is a reduced basis for a segment S if it 
verifies the following properties: 

(i) the polynomials in B are normalized to have content 1 wrt x over A (in 
order to work with polynomials instead of rational functions) ; 

(ii) the leading coefficients of the polynomials in B are different from zero on 
every point of S; 

(iii) B specializes to the reduced Grobner basis of Ca {I), keeping the same lpp 
(leading power product set) for each a e S, i.e. its lpp set remains stable 
under specializations within S. 

Reduced bases are not faithful, i.e. they do not, in general, belong to /. They are not 
unique for a given segment, but the number of polynomials as well as the lpp are unique. 

Definition 3 (Reduced CGS). A CGS is said to be reduced if its segments have reduced 
bases. 

As it is known, the lpp of the reduced Grobner basis of an ideal determine the cardinal 
or dimension of the solution set over an algebraically closed field. This is the reason why 
disjoint reduced CGS are very useful for applications as they characterize the different kind 
of solutions of V(I). 

Using Weispfenning's suggestions the author [Mo02] obtained an efficient algorithm 
(DISPGB) for Discussing Parametric Grobner Bases to compute a disjoint reduced CGS. 
Actually this algorithm is called BUILDTREE. 

BUILDTREE builds up a dichotomic binary tree, whose branches at each vertex cor- 
respond to the annihilation or not of a polynomial in K\a]. It places at each vertex v a 
specification = (N V ,W V ) of the included specializations, that summarizes the null and 
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non-null decisions taken before reaching v, and a specialized basis B v of cr-s(I) for the spe- 
cializations cTq £ E„. The set of terminal vertices form a disjoint reduced CGS where the 
segments S v are characterized by reduced specifications determined by (N V ,W V ). 

Since then, more advances have been made. Inspired by BUILDTREE, Weispfen- 
ning [We03] gave a constructive method for obtaining a canonical CGB (CCGB) for para- 
metric polynomial ideals. Using this idea, Manubens and Montes [MaMo06] improved 
BUILDTREE showing that the tree To built up by BUILDTREE, can be rewritten as a 
new tree T providing a more compact and effective discussion by computing a discriminant 
ideal that is easy to compute from To. The rebuilding algorithm given in [MaMo06] can 
be iterated to obtain a very compact new tree organized as a right-comb tree. It builds an 
ascending chain of discriminant ideals that orders the segments defined as differences of the 
varieties of two consecutive discriminants, and provides a very compact disjoint reduced 
CGS. Nevertheless this rebuilding algorithm does not always produce the canonical CGS. 
This method will be presented in a forthcoming paper. 

Having in mind the improvement of our disjoint reduced CGS to obtain a canonical 
CGS, in the present paper we adopt a different perspective. Instead of rebuilding the tree, 
we analyze all the different situations that can occur in the BUILDTREE CGS, formulate 
a natural conjecture and, using it, show how the segments can be packed to obtain the 
largest possible segments allowing the same reduced basis. These segments become non- 
algorithmic dependent, i.e. intrinsic for the given ideal. The packed intrinsic partition 
contains the minimum number of segments corresponding to reduced bases. Algorithms to 
perform the discussion about which segments must be packed, to obtain the reduced basis 
for the packed segments and to describe segments defined by difference of two varieties in 
a canonical form are also given. It remains to give a canonical description of the union of 
the segments included in the packed segments as well as the algorithm to carry it out. This 
last step is described in [MaMo07a]. 

The whole set of algorithms, including those in |MaMo07a] have been implemented^ and 
denoted MCCGS (Minimal Canonical Comprehensive Grobner System) algorithm. It is yet 
operative and in experimental phase. It is promising and very useful for applications as 
can be seen through its applications to automatic geometric theorem proving and discov- 
ering |MoRe07j . It may be objected that canonicity pays a price in computing time. The 
implementation shows that in fact the time increases only about 20-30% whether the output 
becomes much more simpler, compact, easy to be understood and practical for applications. 

Recently Sato and Suzuki |SuSa06] have developed a new simple algorithm for computing 
CGS based on Kalkbrenner's Theorem [Ka97j , Its interest lies in its simplicity (it is perhaps 
sometimes more efficient) but the output is not sufficiently clear and useful for applications. 

Section [2] reviews the basic features of BUILDTREE. In Section [3] it is explained what is 
meant by canonical CGS, a conjecture is formulated and using it, it is shown how to obtain 
the canonical CGS. Section 0] give the basic theorems to obtain a canonical description 
of diff-specifications and gives the corresponding algorithms. Finally, in Section [5] further 
developments are pointed out giving some insight as an advance of the content of the final 
paper [MaMo07a], where the whole canonical description of the union of segments will 

Manubens and Montes implementation of the new algorithm MCCGS is available on the web 
http://www-ma2.upc.edu/~montes. The library, called DPGB release 7, is implemented in Maple 8. 
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be given making the MCCGS the algorithm to provide a canonical representation of the 
intrinsic segments. 

In this paper, we only give partial examples that illustrate the algorithms discussed here. 
A unique complete example is given in the final Section [5] to give an idea of how simple is 
the final output that has the minimum number of segments. 

The algorithms work with ideals, as these are the algebraic objects that allow a Grobner 
representation. But ideals do not represent varieties in a unique form. So we frequently 
adopt a geometrical view. We shall consider ideals defined in A = K\a], where K is the 
computable field (for example Q), whereas the varieties will be considered in K' m , where 
K' is an algebraically closed extension (for example C). Let J be an ideal in K\a], and 
J' = J ■ K'\a] be its extension to -RT'fa]. The symbol V(J) will denote 

V(J) = {a G K' m : V/ G J, /(a) = 0} 

= {a G K' m : V/ G J', /(a) = 0} = V(J'). 

We emphasize the use of the non-standard notation V in the whole paper as used in the 
extended affine space, whereas the ideals are defined in the base field. Lemma[8]in SectionH] 
will justify that decision. 



2 Reviewing BUILDTREE 

A specification of specializations is a subset £ of specializations determined by a con- 
structible set of the parameter space. 

BUILDTREE uses reduced specifications for the segments. Different definitions have 
been given in [Mo02] and [MaMo06]. The reduced specification used in release 4 of the 
DPGB package described in [Mo02] does not require N to be radical nor to obtain a prime 
decomposition of N. In this approach, when we need to test whether a polynomial in K\a] 
vanishes for a G E, it is not sufficient to divide it by N. Instead, we must test if it belongs 
to y/ (N). But this is simpler than computing the radical and its prime decomposition. This 
makes REDSPEC more efficient but does not ensure all the nice properties that we want 
to have. Nevertheless, even if this is a good practical solution, for theoretical purposes we 
need to replace the concept of reduced specification^. 

Definition 4 (Red-specification). Given the pair (N,W) of null and not null conditions 
denote 

h= Y[w£ K[a], and Y{h) = [j Y(w) C K' m . 

They determine a reduced specification of specializations (red- specification) whenever 

1. N is a radical ideal described by its reduced Grobner basis wrt 

2. W is a set of distinct irreducible polynomials in K\a], 

3. Let Ni be the prime components of (N) over K\a]. Then h iVj for all i. 



definition 7 in MaMo06 . 
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Note that properties (ii), (iii) of the definition in MaMo06] are simple consequences 
of Definition |U Nevertheless property (3) of Definition 2] is stronger, and REDSPEC (de- 
noted CANSPEC in previous papers) is supposed here to verify this new definition of 
red-specification. 

The segment associated to a red-specification is SWw) = ^(-^0 \^(M> an d the included 
specifications are £(jv,w) = i a a '■ a £ Y(N) \ Y(h) G K"™}. Let W = {u>i, . . . ,uj s } C K\a] 
and A = (Ai, . . . , A s ) G Z> - Define 

s 

W(A)=aJ X = n^ Ai - 

i=l 

The set of all non-null polynomials of K\a] as a consequence of is 

W* = {kW(X) : k G K, A G (Z^ ) s }. 

Definition 5 (Reduced polynomial). A polynomial / is reduced over the segment S deter- 
mined by the red-specification (N, W) if / = /, con%(/) = 1 and lc(/) G W*. 

Definition 6 (Good specialization). We say that the polynomial F specializes well to the 
reduced polynomial / over the segment S determined by the red-specification (JV, W), if F 

— TV 

and / are proportional except for non-null normalization, i.e. if aF = bf with a,b G W*, 
(i.e. the coefficients a, b do not become on any point of S). 

BUILDTREE is a Buchberger-like algorithm. Applied to the ideal / it builds up a 
rooted binary tree with the following properties: 

1. At each vertex v a dichotomic decision is taken about the vanishing or not of some 
polynomial p(a) G K\a\. 

2. Each vertex is labelled by a list of zeroes and ones; the root label is the empty list. 
At the null child vertex p(a) is assumed null and a zero is appended to the parent's 
label, whereas p(a) is assumed non-null at the non-null son vertex, in which a 1 is 
appended to the father's label. 

3. At each vertex v, the tree stores (N v , W v ) and B v , where 

- (N v , W v ) determines a reduced specification T, v of the specializations summariz- 
ing all the decisions taken in the preceding vertices starting from the root. 

- B v is reduced wrt Y> v (not faithful) and specializes to a basis of a-a (I) for every 
^ao ^ S 0! preserving the 1pp. 

4. The set of terminal vertices form a disjoint reduced CGS in the sense of Definitions Q] 
and[3j 

- B v specializes to the reduced Grobner basis of ctq (7) for every Oa G X„ and 
has the same lpp set. The polynomials g in the bases are normalized having 
con%(g) = 1. 
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- The specifications of the set of terminal vertices ti determine subsets C K' m 
forming a partition of the whole parameter space K' m : 

X = \St i ^iu • • • > St p }, 

and the sets have characteristic lpp sets that do not depend on the algorithm. 

5. The unique vertex having as label a list of 1 ([1, . . . , 1]) corresponds to the generic 
case as it is determined by only non-null conditions. It does not necessarily contain 
the whole generic case, as we will see next. 

Thus the terminal vertices of BUILDTREE form a disjoint reduced CGS. 

3 Finding a canonical CGS 

The objective of the paper is to advance in the definition and computation of a unique 
(canonical) CGS. In order to reach this objective we need to obtain an intrinsic family of 
subsets Si for the CGS, uniquely determined. 

Denote T = (d,...,C s ) the disjoint reduced CGS built by BUILDTREE, i.e. the list 
of terminal cases C% = (B{, Si). We shall always set the generic case as the first element of 
T. We group them by lpp. 

r = ((Cn, . . . , Ci si ), . . . , (Cfci, . . . , Cks k )) = (ri, . . . , r^) 

where the first index denotes the lpp and, as usual, C\\ corresponds to the fundamental 
segment of the generic case. Obviously the sets in each group (Cu, . . . , Cj Si ) are canonically 
separated because it cannot exist a common reduced basis for them, as reduction implies 
preservation of the lpp. Thus if it is possible to obtain a unique reduced basis for each 
group then we will have a canonical CGS. But even when this is not possible and one 
or more groups must be split into several subgroups forming canonical equivalence classes 
where each class admits a common reduced basis, we will have a canonical CGS. Thus our 
objective is to obtain this classification. 

Conjecture 7. Let Ti = (Cu, . . . , Ci Si ) be the set of all segments of a disjoint reduced CGS 
having reduced bases with a common lppj . If Cij and Cn~ admit a common reduced basis 
and Cik and Cu also, then it exists a common reduced basis to Cij, Cn~ and Cu. 

If the conjecture is true then we have an equivalent relation between the segments in Tj 
that is independent of the algorithm and thus shows the existence of the canonical CGS. 
This canonical CGS will be also minimal in the sense that it contains the minimum number 
of segments of a disjoint reduced CGS. We are now concerned with the task of giving 
algorithms to carry out the task of summarizing the Cij varying j forming the equivalent 
class with a unique reduced basis. 

3.1 Using sheaves 

Before tackle that task we need to know that in some special cases we will need to use 
sheaves. We are indebted to Michael Wibmer [Wi06] for the idea of using sheaves for 
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not null 



null 




[1] [1] [1] □ 



Figure 1: BUILDTREE for I = (ax + b, cx + d). 



summarizing some kind of segments, as they are needed for some special problems. Let us 
give an example from him. 

Exemple 1. Let I = (ax + b,cx + d). Applying BUILDTREE with = lex(x,y) and 
>~a = lex(o, b, c, d) we obtain the tree of Figure [TJ that provides the following segments 
(ordered by lpp): 



lpp 


basis 


null cond. 


non-null cond 


[1] 


[1] 


[] 


{a, ad — cb, c} 


[1] 


[1] 


[c] 


{a,d} 


[1] 


[1] 


[a] 


{b,c} 


[1] 


[1] 


[c, a] 


{b,d} 


[1] 


[1] 


[d, c, a] 


{b} 


[1] 


[1] 


[c, b, a] 


{d} 


[x] 


[cx + d] 


[ad — cb] 


{a,c} 


[x] 


[ax + b] 


[d,c] 


{a} 


[x] 


[cx + d] 


[b,a] 


{c} 


[} 


{) 


[d, c, b, a] 


{} 



Obviously, the six cases with basis [1] can be summarized into a single case. We must add 
the six corresponding segments, and thus we will need a method to do this in a canonical 
form. But in any case, the union of the six segments is intrinsic to the problem and 
corresponds to the total generic case having basis [1]. You can see in Section [5] how these 
segments are grouped in the canonical tree build by the MCCGS algorithm. 

The three cases with lpp = [x] can also be summarized into a unique basis but now 
instead of a single polynomial we must use a sheaf with two polynomials. Effectively, 
the polynomial cx + d specializes well in the first and third segments with lpp = [x] and 
specializes to in the second segment. And the polynomial ax + b that forms the reduced 
basis of the second segment is proportional (and thus equivalent) to cx + d in the first 
segment, but specializes to in the third segment. The common basis for the three segments 
in this case is given by one sheaf [{cx + d,ax + b}] instead by a polynomial: at least one 
of the two polynomials in the sheaf specializes well in the union of the segments whether 
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the other either specializes also or goes to zero. Thus we see that for our objective we 
must admit sheaves also for the bases instead of single polynomials. The three segments 
are grouped in the canonical tree. 

When a reduced basis of a segment S contains a sheaf, then we need that, for all a G 5, 
at least one of the polynomials in the sheaf specializes to the corresponding polynomial of 
reduced Grobner basis of the specialized ideal and the others either specialize also, either 
to it or to 0. 

Thus the canonical CGS for this example will contain only three segments, namely 



lpp 


basis 


sets of pairs (N, W) 


[1] 


[1] 


([],{«, ad- cb,c}), ([c],{a,d}), ({a],{b,c}), 
([c,a},{b,d}), ([d,c,a],{b}), ([c,b,a],{d}) 


[x] 


[{cx + d,ax + b}] 


([ad - cb],{a,c}), ([d,c],{a}), ([b,a],{c}) 


[} 


[} 


([d,c,b,a},{ }) 



Sheaves will appear only in over-determined systems with generic basis [1]. For these 
kind of systems, and for a combination of the parameter values making compatible the 
redundance with some degree of freedom as is the case in the previous example, sheaves 
may appear. Nevertheless this is not so for other kind of systems. An example having a 
larger sheaf for the basis of one of his segments is I = (ax + 6, cx + d, ex + /) . 

3.2 Obtaining common reduced bases 

In most common situations where the BUILDTREE CGS presents multiple segments with 
the same lpp it will exist one subsegment (the most generic one) whose basis already spe- 
cializes well in the other segments, and then we only have to pack them. 

There are also problems where it does not exist a common basis for segments having 
the same lpp. A new example from Wibmer [WLQ6] shows this situation. 
Exemple 2. Consider the following simple system: / = (u(ux + l),(ux + l)x). RE- 
BUILDTREE gives the following GCS with two unique segments having the same lpp. 



lpp 


basis 


null cond. 


non-null cond 


[x] 


[ux + 1] 


[] 


{«} 


[x] 


[x] 


[u] 


{} 



It is easy to convince oneself that it does not exist a common reduced basis for both 
segments as the leading term of the generic segment specializes to for u = whether the 
independent term is always different from zero. 

A fourth possibility arises when we have two segments with the same lpp sets, charac- 
terized by (Bi,Ni, W\) and (B 2 , N 2 , W 2 ) that do not directly specialize one to the other by 
reducing the basis, but nevertheless it can exist a more generic reduced basis specializing 
to both. Let us explore that possibility. We want to test if it exists a basis B\2 such that 

(i) lpp(B 12 ) = lpp^) = lpp(B 2 ). 
(h) a (NuWl) (B i2 ) = B 1 and a {N2)W2) (B 12 ) = B 2 
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F <— GENIMAGE(/i, N±, W±, /2, N 2 , W2) 
Input : 

(/1 = ^2 a a x a , N±, W\): basis and red-spec of a terminal case 
(/2 = b a x a , N2, W2): basis and red-spec of a terminal case 
L, M G Z bounds for the tests. 
Output : 

F: when it exists, F returns a polynomial such that C(Ni,Wi)(F) = fi 
and cr(Ar 2 ,w/ 2 )(-F) = fi else it returns F = false. 

begin 

test:= false; F := false 

N := GBEX(A r i + N2) (returns the Grobner basis and also the matrix M 
expressing the polynomials in N in terms of the polynomials in N\ and N2) 

for all A G Z| , |A| < L while not test do iwi = Wi(A) 
for all Jl g Z> , < M while test do tt>2 = W2(ji) 
{ HT is the index of the leading term) 

h := kiWidHT - k2W 2 bHT 
test := true 

if h has a factor Afci + Bk2 with A, S G K then 
A^-^ * == .S . • ~~ -^4. 

for all terms a of /1 or /2 while test do 
/i := k[wia a — k 2 W2b a 
r := h 

if r / then test := false end if 
end for 
if test then 

F:=0 

for all indices a of terms of /1 or f2 do 
h := k[wia a — k' 2 W2b a 

—N 

qi := list of quotients of the exact division h 

I \N\ |JVi| 

F := F + k' lWl a a N V<H M H 

\ i=l 3=1 

end for 
else 

test := false 
end if 
end do 
end do 
end 



in order that both cases can be summarized into a single one conserving the 1pp. Denote 
fi G B\ and /2 G B2 two corresponding polynomials with the same 1pp. Then we must test 
if it exists a F\2 such that a^i,Wi )(^12) = fi and cr(jv 2 ,W2)(-^i2) = fi- For this it must exist 
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wi G W*, n\ G (Ni) ■ K[x] and w 2 G W 2 *, n 2 G (iV 2 ) • K[x] such that 

-F12 = u>i/i + ni = ^2/2 + n 2 . 

Let /i = S Q o Q x a , / 2 = X a 6 a x a and F12 = T, a c a x a . For every index a of a term in /1 or f 2 
the coefficients must verify 

c a = w\a a + n\ a = w 2 b a + n 2a . 

with fixed G W* and 102 G W 2 * and appropriate ni Q G (N\) and n2 Q G (-/V2). This implies 

wia a - w 2 b a = n 2a - n\ a G (Nt) + {N 2 ) = (N) 

that can be solved by testing for all the possible w\ G W* and w 2 G W 2 * when it exists a 
left hand side belonging to (N). The semi-algorithm GENIMAGE does it and will obtain 
an F\2 if it exists. It is a semi-algorithm in the sense that the possible choices of w\ and 
w 2 are in fact not finite and the algorithm must set a bound on the possible total degree 
(\M = Ylt — L) of the terms tested for which no bound is known. Even if this depends on 
a luck, and little combinatorics is used, in practice it does not cause big problems because 
this does not occur often and when it does the result is easily found. 

The semi-algorithm is self understanding. Also, when two or more segments given by 
red-specifications {N\,W\), . . . , (N S ,W S ) have been generalized to a generic basis Bq and 
we must test if a new segment (B, (N,W)) admits a common pre-image, we can also use 
GENIMAGE for each polynomial /q G Bq and / G B taking for /o as null and non-null 
common conditions (n|iVj, DfWj). If GENIMAGE obtains a pre-image it will reduce well 
to all the segments. Let us give an example: 

Exemple 3. Consider the following example from Sato-Suzuki [SaSu03] : I = (ax 2 y + a + 
3b 2 ,a(b — c)xy + abx + 5c) wrt >~x = lex(x, y), >~a = lex(a, b, c). REBUILDTREE obtains 
a CGS with three segments with basis [1] that obviously can be directly added, three cases 
with lpp set [y, x] that do not specialize one to the other, and five other segments with 
distinct lpp namely [y 2 ,x], [y,x 2 ], [yx,x 2 ], [yx 2 ], [ ]. 

The question arises for the three segments with [y, x] as lpp. Let us detail these segments: 



lpp 


basis 


null cond. 


non-null cond 


[y,x] 


[y, 36 3 x - 5c] 


[0 + 36^] 


{b - c, c, 6} 


[y,x] 


[a 2 y + 25, 5x + a] 


lb] 


{c,a} 


[y,x] 


[25y + 3ac 2 + a 2 , ax + 5] 


[b-c] 


{c, a} 



We can verify that none of the bases reduces to the others. Applying GENIMAGE first to 
(£1,(^1,^1)) and (B 2 ,{N 2 ,W 2 )) and then to (#12,(^12 = [b(a + 35 2 )], W 12 = {c})) and 
(S3, (N3, W3)) a common reduced basis is found 

B 123 = [(256c - 25a 3 6 - 756 3 a 2 + 25ca 3 + 75a 2 b 2 c)y 

-625a6 + 1875c& 2 + 625ac + a 2 bc + 3a6 3 c - 18756 3 , 
(ba 2 - 15ab + 15ac - 9b 5 + 96 4 c - 456 3 + A5cb 2 )x 
-3ba 2 + 3a 2 c + 5ab + 276 5 - 276 4 c + 156 3 - 156c 2 ] 

that reduces to the three bases in the respective segments. 

We have explored four possible situations for a pair (or a collection) of segments (Bx,Ni, W±) 
and (B 2 , N 2 , W 2 ) with the same lpp, namely 
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1. the polynomials of B\ reduce to the polynomials of Bi on {N-iiW-i), (most frequent 
case); 

2. first case does not happen but it exists a pre-image basis B\2 that reduces to both 
and can be computed by GENIMAGE; 

3. both bases can be summarized using sheaves; 

4. a reduced common basis does not exist. 

Table [1] shows the algorithm DECIDE that decides if two corresponding polynomials of B\ 
and B2 have a common pre-image or a sheaf or it does not exist. If S(f\, fa) 7^ then it 

calls GENIMAGE that will decide if a pre-image exists or not, but in this case the result 
N2 

cannot be a sheaf. If S(fi,fa) = as lc(/i)/2 — lc(/2)/i specializes to in the subset 

N 2 

S2, then /1 specializes either to fa or to in S2. Then, if lc(/i) G W 2 * then it is always 
non-null in S2 and so f\ specializes to fa and is the generic polynomial F. Else we carry out 
the symmetric comparisons and conclusions, and only when the S-polynomial specializes 
to both in S\ and in £2 and none of the leading coefficients remain non-null in the other 
segment we will have a sheaf {ft, fa}. 

4 Canonical specifications 

The following Lemma plays an important role in the obtention of canonical specifications 
of diff-specifications. 

Lemma 8. Let K be a field of characteristic zero and K' an algebraically closed extension, 
P and Q ideals in K\a], P prime and Q <£_ P. Then, on K' m 

V(P)\V(Q) =V(P). 

Proof. Denote P' = P ■ K'\a] and Q' = Q • K'\a] the respective extensions of P and Q in 
if'[a]. To prove the lemma we follow four steps: 

(i) As P is prime and Q <f_ P, we conclude that P : Q = P. We 
leave the proof as an exercise. 

(ii) (P : Q)' = P' : Q' . See [ZaSa79| . Vol II, p. 221. 

(iii) As P is prime, P' is radical. See [ZaSa79j .Vol II, p. 226. 

(iv) Since K' is algebraically closed and P' is radical, 

Y(P') \ V(Q') = V(P' : Q'). 
See |UoLiSh92j . Theorem 7, p. 192. 

Combining these four steps, we obtain 

V(P) \ V(Q) = V(P') \ V(Q') = V(P' : Qf) = V((P : Q)') = V(P') = V(P). 

□ 

Using Definition [J] we can now prove the following 
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F <- DECIDE(/i ,Ni,W 1 ,f 2 ,N 2 ,W 2 ) 
Input : 

(fi,Ni, W\): basis and red-spec of a terminal case 
(f 2 ,N 2 , W 2 ): basis and red-spec of a terminal case 
Output : 

F: if it exists a pre-image or a sheaf then F is a polynomial (or sheaf) such that 

u {Nx,Wx){F) = fi and (J (N 2 ,W2)( F ) = h else it returns false 
begin 

if5(/i,/ 2 ) JV V0then 

F := GENIMAGE(/i,iV 1 ,M/i,/2,A r 2, W 2 ) 
else 

JVb 

if lc(/i) 6 W 2 * then F := /i 
else 

if SU\Jij Nl ^0 then 

F := GENIMAGE(/i,JVi,Wi,/2,JV2,W2) 
else 

if Hf2) Nl e W{ then F := f 2 
elseF:={h,f 2 } # (sheaf) 
end if 
end if 
end if 
end if 
end 



Table 1: 

Theorem 9. Let (N, W) determine a red-specification. Then we have 



Y(N) \ [j V(w) = Y(N) \ Y(h) = V(JV) 

Proof. Decompose sj (N) = f]j iVj into primes in K\a), so that 

Y(N) \ Y(h) = (\j V(Nj\ \ Y(h) = [J (V(JV,) \ ¥(/>)) . 

As (N, W) determines a red-specification, h iVj for all z, and thus, applying Lemma [8] for 
each i it results 

Y(N)\Y(h) = \JWQ\W) = \JV(Ni) = Y(N). 

i i 

□ 
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We have seen that if Conjecture [7] is true it exists an intrinsic canonical partition of 
the parameter space K' m and a reduced basis for each segment and from the BUILDTREE 
output the algorithms DECIDE and GENIMAGE will obtain it. It is apparent the need 
of giving a canonical representation of the intrinsic partition because otherwise we cannot 
verify its uniqueness, for example if determined by another algorithm. So we focus now in 
the canonical description of the union of red-specifications. 

Definition 10 (Diff-specification). Given two ideals A C M whose associated varieties 
verify V(A) D V(M), they define a diff-specification (A, M) describing the subset S = 
V(iV)\V(M) of K' m . 

In particular a red-specification (A, W) is easily transformed into a diff-specification. 
Take h = Y\ weW w and M = A + (h). Obviously (A, M) is a diff-specification. 

We begin giving a canonical representation of the subsets of a diff-specification, and 
then we shall discuss how to add subsets defined by diff-specifications. 

Definition 11 (Can-specification). A can-specification of a subset C is a representation 
defined by the set of prime ideals (Aj, My) varying i,j such that 

C = Y(N) \ V(M) = [J (V(A t ) \ (U^My )) , (1) 

i 

where N = n,Aj and = HjMjj are the prime decompositions over K\a] of the radical 
ideals M and M.i respectively, where A, C A4ij. 

We have the following 

Theorem 12. 

Every set C = Y(N) \ Y(M) C A' m corresponding to a diff-specification admits a 
can-specification, and the algorithm DIFFTOCANSPEC given in Table \M builds it. 

2. Over K' m , a can- specification verifies 

C = \J (Y(Ni) \ (U,-V(My)) = |J V(iVi) = V(A0. 

i i 

3. The can-specification associated to a set C given by a diff-specification is unique. 

4. All points inCD V(JV f ) are in V(A;) \ (U j V(M ii )). 

Proof. 1. Let y/N = Hi-^i De ^he prime decomposition of the radical ideal vJV over 
K\a}. Then we have 

C = Y(N) \ Y(M) = (\J V(JV<) J \ V(M) = (J (V(A,) \ V(M + A,)) . 

In this decomposition the variety to be subtracted from V(Aj) is contained in it. 
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S <- PRIMEDECOMP(iV) 
Input : 

N: ideal (representing a variety) 
Output : 

S = (Ni, ... ,Nk): the set of irredundant prime ideals wrt K of the decomposition 
oi^N = HjNj 

Y «- DIFFTOCANSPEC(iV, Af) 
Input : 

N: the null-condition ideal of the diff-specification 
M: the non-null condition ideal of the diff-specification M N 
Output : 

Y = {(Ni, ({Mij : 1 < j < li})) : 1 < i < k}: the set of prime ideals corresponding 
to the canonical decomposition of Y(N) \ V(M) (Theorem [T2|) 

begin 

Y = 

S := PRIMEDECOMP(iV) 
for Nj G S do 
if iVj / y/M + Nj then 
Tj- := PRIMEDECOMP(M + Nj) 
Y^YUjiiNj,^)} 
end if 
end for 
end 



Table 2: 

It can happen that y/M + Ni = (1), in which case nothing is to be subtracted from 
Y(Ni). It can also happen that \JM + N] = N { , in which case the term V(A r i ) \ V(JVj) 
disappears from the union. The above expression is simplified and for all the remaining 
terms we have Ni C y/M + N { 

Let now \/M + JV ; = flj M i:; - be the prime decomposition of y/M + Ni over K[a]. For 
each j we have Ni C M«. The decomposition becomes 

V(N) \ V(M) = J (V(JVi) \ (Uj-VCAfy)) , (2) 

i 

where N = f] t N and Aii = P| ■ M^- are the prime decompositions over K[a] of the 
radical ideals M and M.i respectively, proving part (i) of the theorem. (Observe that 
the algorithm DIFFTOCANSPEC is nothing else than the description done in this 
paragraph). 

It should be noted that the prime decompositions in the computations are performed 
in K\a], as it is the computable field. Thus these decompositions can split over K'\a]. 
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In the same sense we cannot ensure that the varieties V(iVj) nor V(M,j) are irreducible 
over K m nor over K' m as we cannot use the Nullstellensatz in K. Nevertheless the 
prime decompositions are canonically well defined over K\a}. 

2. Using Lemma[8]for each term in the decomposition given by formula ([I]) of C we have 

C = {J V(iVi) \ (UfV(My) = |J Y(Ni) = Y(N) 

i i 

over K' m , proving part (ii) of the theorem. 

3. Suppose that C admits two diff-specifications characterized by the pairs of discrim- 
inant ideals (N, M) and (R, S) respectively. If we denote by 1Z = f\ Rk and Sg = 

Ski the ideals in the decomposition obtained from the diff-specification with R and 
S using the method described in part (i) of this theorem they will verify C = Y{J\f) = 
V(<S) by part (ii). As J\f and S are radical, they are also radical over K'[a] and thus 
by the Nullstellensatz they are both equal to 1(C) over if [a]. Thus we have N = Ri 
for each i, as the prime decomposition in K\a] is unique. 

Next we subtract from each V(iVj) all the points that are not in C as they are in Y{M). 
We have already eliminated the components of Y(N) that are also in Y(M), so that 
the points in V(iVj) that are not in C are the points of the variety V(M + iVj) C V(iVj). 
Then by the Nullstellensatz Mi = y/M + Ni is the variety ideal I(Y(Mi)). Carrying 
out the prime decomposition of Mi we are done with the canonical decomposition. 
Thus the decomposition of C given in part (i) of formula ([T|) is unique. 

4. This is now obvious as we have subtracted from each Y(N) all the points in Y(N) n 
V(M). 

□ 

Note that the can-specification is canonical but the constructible sets whose union de- 
scribes C do not have empty intersection. Nevertheless this does not affect C itself. 

5 Further developments 

If Conjecture [7| is true, it exists a minimal canonical CGS. The algorithms here described, 
start from the BUILDTREE CGS and regroup the segments to obtain the intrinsic segments 
having a reduced basis. The result for each segment is of the form: 

d = (B h Si) = (Bi, ((N a , W a ), (Ni^WijJ)) 

The subsegments defined by (Nik,Wik) are described by red-specifications and thus as a 
difference of varieties V(iV^) \ V(M^) where h = Y\ weW . w and = {N^) + (h) thus 
{Nik, Mik) is its diff-specification. In order to obtain a canonical description of the intrinsic 
partition of our disjoint reduced CGS it is apparent that we need to add diff-specified sets 
in a canonical form. This task is done in [MaMo07a] . 
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[d, cjb, a] 



Figure 2: MCCGS tree for Exemple 1. 



Let us outline how this works. We cannot assume that the simple form given by for- 
mula (pQ) will be sufficient. A more complex constructible set will be formed. There can 
exist different canonical forms for describing it, but in any case this will need prime decom- 
position of radical ideals. Our canonical form is given by an even level rooted tree called 
P-tree whose root defines level 0. At the nodes there are prime ideals Pi 1) ... ) i j of K\a\. These 
ideals verify Pi lt ... t i- C Pj^...^.^ for every k and the set of Pi lt ... t i ,fc f° r every k are the prime 
decomposition of a radical ideal T 7 ^,...,^. The set C defined by the P-tree has to be read 

c = {jY(p tl )\ (ljv(p ilia )\ (iJ^^WsA (••■ U 7f/ '— j 

Let us advance some results from [MaMo07a] and illustrate how is the final output of the 
MCCGS algorithm for Exemple 1. It is illustrated by a plot procedure in Figure [21 and the 
algebraic output summarized in the following table: 



lpp 


basis 


segment 


[1] 


[1] 


C 4 \ (V(ad - be) \ (V(o, c) \ V(a, b, c, d))) 


[x] 


[{ex + d,ax + b}] 


V(ad - be) \ V(a, c) 


[} 


[} 


Y(a,b,c, d) 



It can be seen that the MCCGS algorithm gives a very compact solution for the problem 
easy to interpret. This is generally so for many other problems. It has been successfully 
applied to geometrical theorem discovery [MoRe07] obtaining very simple answers for rela- 
tively complex problems. 

There are two possible lacks coming from the Conjecture and the semi-algorithm GEN- 
IMAGE as we must set artificial bounds to make it algorithmic. Nevertheless, the use of 
MCCGS is at least useful to find examples where the minimal canonical CGS either does 
not exist or is not obtained by the actual algorithm, providing examples to test both the 
Conjecture and the semi- algorithm. 

Finally it must be pointed out that the term order >~a chosen for the computation in 
A only affect the description by Grdbner bases of the varieties describing the segments but 
not to the varieties themselves. 
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